BPS Solutions to a Generalized Maxwell-Higgs Model 
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We look for topological BPS solutions of an Abelian-Maxwell-Higgs theory endowed by non- 
standard kinetic terms to both gauge and scalar fields. Here, the non-usual dynamics are controlled 
by two positive functions, Gdcjil) and which are related to the self-dual scalar potential 

y(|</)|) of the model by a fundamental constraint. The numerical results we found present inter- 
esting new features, and contribute to the development of the recent issue concerning the study of 
generalized models and their applications. 
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I. INTRODUCTION 

In the context of classical field theories, topological 
structures, such as kinks vortices and magnetic 
monopoles Q, are described as static solutions to some 
nonlinear models. In particular, such models must al- 
low for the spontaneous symmetry breaking mechanism, 
since topological solutions are formed during a symme- 
try breaking phase transition. In this sense, topologi- 
cally non-trivial configurations are of great interesting 
to physics [4] , mainly concerning the cosmological conse- 
quences they may engender, since such configurations can 
appear in a rather natural way, during phase transitions 
in the early universe. 

In particular, vortices are described as rotationally 
symmetric solutions of a planar Abelian-Maxwell-Higgs 
model endowed by a fourth-order scalar potential for 
the matter self-interaction, which also introduces sym- 
metry breaking and nonlinearity. In the usual case, the 
Maxwell term controls the dynamics of the gauge field, 
and the covariant derivative squared term controls the 
dynamics of the scalar field. In such context, vortices 
are finite-energy solutions of a set of two coupled first- 
order differential equations, named Bogomol'nyi-Prasad- 
Sommerfield (BPS) equations In this case, the BPS 
vortices are the minimal-energy solutions of the model, 
and they have interesting applications, mainly concerning 
the superconductivity phenomena and the superfiuid 
He*. 

During the last years, beyond the standard configura- 
tions previously cited, modified ones, also named topo- 
logical k-solutions, have been intensively studied. These 
solutions arise in a special class of theoretical field mod- 
els, generically named k-theories, which are endowed by 
non-usual kinetic terms. As expected, such terms change 
the dynamics of the overall model under investigation. 
Moreover, it is important to reinforce that the idea of a 
non-standard dynamics arises in a rather natural way, in 
the context of string theories. 



In general, k-thcories have been used as effective mod- 
els mainly in the cosmological scenario. Here, many au- 
thors have used the so-called k-essence models Q to in- 
vestigate the present accelerated inflationary phase of the 
universe Also, such models can be used to study 
strong gravitational waves dark matter tachyon 
matter |10| and others [ll| . There are also important mo- 
tivations concerning the strong interaction physics; see, 
for instance, Ref. |12| . 

In this context, the overall conclusion is that the exis- 
tence of topological structures is quite sensible to the use 
of non-standard kinetic terms. Here, the more interest- 
ing issue is that such structures can exist even as solu- 
tions of some theories which do not allow for the sponta- 
neous symmetry breaking mechanism [l^. On the other 
hand, a rather natural way to study such topological so- 
lutions is comparing them with their canonical counter- 
parts. In this sense, some of us have already investigated 
topological solutions in the context of k-field models en- 
dowed by spontaneous symmetry breaking potentials for 
the scalar-matter self-interaction, and interesting results 
can be found, for instance, in Ref. Other interest- 

ing results concerning k-field models and their classical 
solutions can be found, for instance, in Ref. [isj . 

Another interesting issue concerning topological k- 
solutions is that they can be either much larger or much 
smaller than their usual counterparts. In this sense, k- 
bosons can mediate either large-range or small-range in- 
teractions. Also, important physical quantities, such as 
energy density and electric and magnetic fields, can ex- 
hibit proeminent variations on their tipical profiles, in- 
cluding on their maximum values; for a detailed treat- 
ment of some of such features, see Ref. (iGj . 

In the present context, some k-theories can also sup- 
port topological solutions with a finite wavelength. These 
solutions, generically named compactons [17| . are quite 
different from the standard topological structures, which 
interact even if separated by an infinite displacement, 
since they have an infinite wavelength: two adjacent com- 
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pactons will interact only if they come into close con- 
tact, due to their already explained finite wavelength. 
In this sense, compactons are most appropriated to de- 
scribe particle-like configurations than the usual non- 
compactified topological structures. 

So, in this work we present new results concerning the 
topological solutions of an Abelian-Maxwell-Higgs model 
endowed by non-usual kinetic terms to both gauge and 
scalar fields. Here, the non-standard dynamics are in- 
troduced by two positive functions, G(|(/>|) and wd^j), 
which couple with the Maxwell term and with the co- 
variant derivative squared term, respectively. The Euler- 
Lagrange equations of motion of such theory are hard 
to solve, then we focus our attention only on the finite- 
energy solutions of the BPS equations of the model. 
These equations can be obtained by minimizing the en- 
ergy functional of the model, which can be done via 
an important constraint between G(|(/)|), and the 

scalar potential (101) for the matter self-interaction; see 
eqs. p3|) and below. In this case, rotationally sym- 
metric BPS vortices are described as the minimal-energy 
solutions of the cited model, and they engender interest- 
ing features, as explained below. 

This paper is outlined as follows: in the next SecHH we 
introduce the model and develop the theoretical frame- 
work which allows us to get to its Bogomol'nyi-Prasad- 
Sommerfield equations. In Sec lIIIl we prove the con- 
sistence of the theoretical framework previously devel- 
oped by using it to investigate the existence of new BPS 
states. Here, we note that such states are constrained 
to the choices made for G(|(?i|) and w(|(/)|): for any ac- 
ceptable pair of such functions, there is a corresponding 
BPS configuration. Also in Sec lIIIl we show how to use 
the theoretical framework presented in this work to re- 
cover the recover the BPS results concerning the standard 
Maxwell-Higgs model. In Section ITVl we perform the nu- 
merical analysis concerning the new BPS states previ- 
ously presented, we depict the corresponding minimal- 
energy modified solutions and comment on their main 
features. Finally, in Sec|Vl we present our conclusions 
and perspectives. 

From now on, we use standard conventions, including 
natural units system, and a plus-minus signature for the 
planar Minkowski metric: diag(?7'^'') — {-\ ). 



netic matter term. Both G(|(/)|) and are dimen- 

sionless functions to be given below, as functions of the 
amplitude of the scalar field. 

In the case of vortices, it is convenient to deal with di- 
mcnsionless variables. So, for simplicity, we introduce the 
mass scale M, and use it to implement the scale trans- 
formations: ^ x'^/M, (f) M^/^cj), At" M^^'^Af, 
e — > Af^/^e and v — M^^^v, where v stands for the sym- 
metry breaking parameter of the model. In this case, 
we get £g — ^ M^£, with C being the dimensionless 
Lagrange density to be used from now on, which has 
the same functional form of Cg- Moreover, we omit 
the coupling constant related to the scalar matter self- 
interaction. Also, we take e and v as real and positive 
parameters. 

To search for vortex solutions, the canonical proce- 
dure is to deal with the Euler-Lagrange equations of the 
model. In the present case, these equations are 



J' 



and 



wd^d'' 101 + 101 df^w - e^A^A'' 101 w 
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Here, we take = Ff.^F^"" and 

= -2e'^w\cf)\^ A^" 



(2) 



(3) 



(4) 



as the modified 4-currcnt vector. 

The Gauss law for time- independent fields can be writ- 
ten as 

Gdkd^A° + dkA°d^G = -2e^w |0|^ A° , (5) 

where E = We note that the Eq. ([5]) is trivially 

verified by A'^ = 0. So, we fix this gauge and use it from 
now on. 

We look for vortex solutions of the form 



(r, 9) = vg (r) e* 



(a{r) 



(6) 



(7) 



II. THE MODEL 

In this section, we introduce the model. It is described 
by the (2 -I- l)-dimensional Lagrange density 

CG = -\G{\4>\)F,,F^'^ + w{\<t>\)\D,^\'-V{\4>\) . (1) 



Here, r and 9 are polar coordinates, and n = 
±1,±2,±3,... is the winding number (vorticity) of the 
solution. In terms of ^ and ([7]), the Euler-Lagrange 
equations ^ and ([3]) can be rewritten as 



G 



dr"^ 



dG 
dr 



G\ da 
dr 



o 2 2 2 
= 2e V g aw , 



(8) 



Here, F, 



di,An is the usual Faraday field 



strength tensor, D^cj) — 9^ 



iAiid) is the covari- 



ant derivative and y (|0|) is the spontaneous symmetry 
breaking potential. Also, G(|0|) is the "dieletric func- 
tion", and w (|0|) |-D^i0| stands for the non-standard ki- 




1 day dG 
er dr J dg 



(9) 
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Equations ([8]) and ^ are the Euler-Lagrange equations 
of motion to the profile functions a (r) and g (r), respec- 
tively. 

To solve ([8]) and ([9]), we need to specify the model. 
In general, we can do it by choosing non-trivial forms to 
G (g) and w{g). In this case, we have to keep in mind 
that both these functions must be positive, in order to 
avoid problems with the energy of the model; see the 
expression for the energy density (|15p below. Also, we 
need to choose a Higgs potential V [g) which allows for 
the spontaneous symmetry breaking mechanism. 

Before that, we note that the limit w (g) — > 1 leads us 
back to the model studied in p^] , which is supported by 
applications concerning the interaction between quarks 
and gluons In this case, the limit G (g) — 1 leads us 
back to the usual Maxwell-Higgs theory. In this context, 
if we choose the symmetry breaking Higgs potential 



Vs (g) 



2 4 

e t; 



(10) 



the Euler-Lagrange equations and (0) can then be 
rewritten as 



(Pa 1 da 
dr^ r dr 



r, 2 2 2 

2e V g a 



(11) 



of such states is still possible, and it is closely related to 
an important constraint between G (g), w (g) and V {g): 



— y/GV V2ev^wg 



(16) 



As we clarify below, for any positive choices to G (g) and 
w {g) , there is a corresponding symmetry breaking Higgs 
potential V (g) as a solution of In this context, gen- 
eralized first order equations can be found by minimizing 
the energy functional (|T5|l . 

To search for BPS vortex solutions in the modified 
model, we need to know how the functions g (r) and a (r) 
behave, near the origin and asymptotically. Near the ori- 
gin, these functions must avoid singular fields. So, they 
have to behave according to 



g{r 0) and a{r 0) n — 1 



(17) 



where we have fixed n = I, for simplicity. Also, the sym- 
metry breaking vortex configurations must have a finite 
total energy. So, as a condition to make the energy finite, 
the energy density ([15]) must vanish. Then, asymptoti- 
cally, g (r) and a (r) have to obey 



g(r ^ oo) — ^ 1 and a (r — > oo) -> 



(18) 



(Pg ^ Idg 
dr"^ r dr 



a^g 



2 2(2 

e V gig 



1) 



(12) 



According to our conventions, eqs.(fTT|) and (fT2|) are com- 
pletely solvable by the first order differential equations 



dg_ ^ ^ga_ 
dr r 



(13) 



In the next Sec. HI, we use the theoretical framework 
developed in this Section to investigate the existence of 
BPS states in the generalized model. Also, we present 
and comment the resulting numerical solutions. Finally, 
we show how to map the standard Maxwell-Higgs and 
Chern-Simons-Higgs first order equations; even in these 
cases, Eq. ([T6| leads us to a physically different rotation- 
ally symmetric solutions. 



-^=±eV(g2-l) . (14) 
r dr 

The solutions of eqs.(fT3|) and ([T4l) are the Bogomol'nyi- 
Prasad-Sommerfield (BPS) states of the standard 
Maxwell-Higgs model PH)) . These solutions are the well- 
known Abrikosov-Nielsen-Olesen (ANO) ones which 
solve the equations of motion (jlll) and (jl2p by minimaz- 
ing the energy of the resulting vortex configurations. 

In general, for non-trivial choices to G (g) and w(g), 
the equations of motion ([8]) and ([9]) will be much more 
sofisticated than the eas. dTTI) and ([T^ . So, in order to 
get an useful insight about the non-trivial case, we con- 
sider the expression for the energy density of the modified 
vortex solutions: 

From Eq. (jl5p . we note that the presence of non-trivial 
G (g) and w (g) makes it hard to obtain the BPS states 
of the modified model. Even in this case, the existence 



III. NEW BPS STATES 

We now pay due attention to the modified BPS states 
themselves. Here, we choose positive functional forms to 
G (g) and w (g). Then, we solve the constraint (IT6|) to get 
to the consistent Higgs potential V (g), which must allow 
for the spontaneous symmetry breaking mechanism. A 
posteriori, we use these conventions to obtain the first or- 
der differential equations of the modified model, by min- 
imizing its energy functional (jlSp . Then, in the next 
Sec. IV, we numerically solve these equations, according 
to the boundary conditions (IT71) and ^TE\\ . Finally, in 
Sec.V, we comment the main features of the resulting 
solutions. 

The model to be studied here has two unusual func- 
tions, G {g) and w (g), which we have to fix to determine 
how they affect the vortex solutions. The function G (g) 
stands for a "dieletric function", and it controls the non- 
standard kinetic term to the gauge field. On the other 
hand, the function w (g) controls the non-usual kinetic 
scalar matter term. Both these functions are dimension- 
less, and are functions of the amplitude of the scalar field. 
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Also, they must be positive, in order to avoid problems 
with the energy of the model; see Eq. (ITSt . 

A rather natural way to investigate the consistence of 
the theoretical framework developed in the previous Sec. 
II is to study the BPS states of the standard Maxwell- 
Higgs model, which is defined by G (g) = w (g) = 1. In 
this case, Eq. leads us to the potential (ITU|) . Then, 
the energy density pSI) can be rewritten as 



1 / 1 da 



da , d , y . 
T—-r±—-r{9 a) . 
r dr r dr 



do qa 
dr r 



(19) 



The resulting total energy E is minimized by the first 
order equations (flBl) and p4)). In this context, the energy 
of the BPS states is 



e (r) d^r = 2ttv'^ \n\ 



(20) 



As usual, this is the lower bound of the energy func- 
tional, i.e., the Bogomol'nyi bound. Also, it is quantized 
according to the winding number n. 

Now, we introduce an interesting new model, which is 
defined by 



G{g) 



(.9^ + 3) 

9 

9 



and w (g) = 2 (g^ + l) . (21) 



According to these choices, we solve (ITCl) to get to the 
potential 



V (g) = g'Vs (g) , 



(22) 



which allows for the spontaneous symmetry breaking 
mechanism, as desired. We point out that this new model 
is not a parametrization of the standard Maxwell-Higgs 
one, since the vacuum manifold of the two models are 
quite different: it is a circle for the usual model ([T0|. 
while it is a dot surrounded by a circle for the new model 



Using eas. ipTI) and (f^ . the energy density ^TE\i can 
then be rewritten as 

_ {g' + 3)' (Ida evYil^y 
^ 2^2 [erdr^ {g^ + 3) / 



dg_^ga 
dr r 



d 



±vi(«(5'-l) (.9' + 3)) 



(23) 



Then, the corresponding energy functional is minimized 
by the equations 



d£ _ _^ga 
dr r 



1 da , e^v'^g'^ (.g^ 



r dr ~'~ ((/^ + 3) 



(24) 



(25) 



We note that the modified equations (p4)) and (|25p can 
not be parametrized into the standard ones, i.e., eqs. (|13|) 
and (fT4| . Even in this case, the energy of the modi- 
fied field solutions is bounded from below, and the Bo- 
gomol'nyi bound is 



(26) 



Here, as in the standard case, the total energy of the BPS 
states is quantized; see Eq. (PH)) . 

In order to reinfoce the consistence of the theoretical 
framework presented in this work, we introduce another 
modified model. It is defined by 

G{g) = {g^ + iy and w {g) = 2g' . (27) 

According to these conventions, Eq. leads us to the 
standard Higgs potential (ITO|) . So, the modified model 
(P7| has the same vacuum structure as the usual model. 
Even in this context, the new model is not a parametriza- 
tion of the standard one; see the BPS eqs. (P^ and ([50)) 
below. The modified energy density can be written in 
the form 




dr 



d 



(28) 



and the resulting energy functional is minimized by the 
first order equations 



d9_ ^ j^ga 
dr r 



1 da 

r dr 



±- 



ia' - 1) 



(5^ + 1) 



(29) 



(30) 



Then, the Bogolmol'nyi bound is also given by (|20p , and 
we conclude that the modified solutions ([^ and ([50)) 
have the same energy of the standard ^3]) and ([M]) ones. 

To end this Section, we show how to use the theoretical 
framework developed in this work to map the standard 
Maxwell- and Chern-Simons-Higgs first order differential 
equations. We do it by reviewing the model studied in 
p^ . which is defined hy w{g) = 1. In this case, as pre- 
sented in that work, an interesting choice to the dieletric 
function is 



(31) 



G(5) = (l + A) 1 + 2A 



Here, A is a real auxiliary parameter which controls the 
model, and k stands for the coupling constant corre- 
sponding to the Chern-Simons term. In the limit A — 0, 
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Eq. ([3T|) gives G{g) — )• 1, and ([16]) reproduces the stan- 
dard Maxwell- Higgs system ([TU]). On the other hand, the 
limit A — oo gives 



G[9) 



We then solve Eq. to get 



y{9) 



-9^9' -ly 



(32) 



(33) 



which is exactly the sixth-order symmetry breaking Higgs 
potential one finds in the usual self-dual Chern-Simons- 
Higgs model. Now, from eqs.(|32|) and ([33|), the energy 
density (llSp can be written in the form 



dr r 



1 da 2e3u4 
er dr 



9 [9 

± w -— (a I 
r dr 



1)) . (34) 



The resulting total energy is minimized by the differential 
equations 



1 da 

r dr 



dg ^ ^ga 
dr r 



fc2 



(35) 



(36) 



and the Bogomol'nyi bound is E = Eg. The equations 
(1551) and ([55]) are just the BPS ones for the standard 
Chern-Simons-Higgs system In this sense, Eq. (pij) 
works as a unified way to map both the Maxwell- and 
the Chern-Simons-Higgs first order differential equations. 
Finally, we point out that, despite the possibility of re- 
producing such Chern-Simons equations, the unified pic- 
ture (|3T]) leads us to a physically different situation, since 
the modified model we consider here only supports non- 
charged field solutions; see the Gauss law ([S]) above. In 
this case, the modified nontopological field solutions are 
unstable, and they can decay into the elementary mesons 
of the model. This fact makes an important difference, 
since the nontopological configurations presented by the 
self-dual Chern-Simons-Higgs theory are stable, and can 
not decay. 

In the next Section, we present the modified first order 
numerical solutions for the profile funtions g (r) and a (r) . 
Also, we comment on the main features of such solutions. 



IV. NUMERICAL SOLUTIONS 

Let us focus attention on the modified numerical so- 
lutions themselves. The equations to be solved are the 
first order ones jMl) and and ^ and Also, 

we solve the standard equations and (|14p . for com- 
parison. In all cases, the functions g (r) and a (r) must 
obey the boundary conditions P7|) and p^ . 



In the present work, the numerical strategy to be em- 
ployed is the relaxation one. In this case, it is necessary 
to input an approximated field solution. Then, our algo- 
rithm will " relax" it into the correct one. To start the nu- 
merical analysis, we consider a variation of the standard 
Maxwell-Higgs model, in which the potential for the self- 
interaction for the scalar field is not present. Then, we 
use its solutions to solve the self-dual Maxwell-Higgs case, 
i.e., the equations (fT3]) and p^ . from which we get the 
well-understood Abrikosov-Nielsen-Olesen (ANO) con- 
figurations. We then use these solutions to initialize the 
numerical study of the modified theory. 

Starting from such ANO configurations, we numeri- 
cally solve the modified BPS equations presented in the 
previous Section, i.e., eas. ([M)) and (pSj). and also eas. ip^ 
and (|30p . for e = u = n = 1. The solutions for the profile 
functions g (r) and a (r) are plotted in Figs. 1 and 2, 
respectively. Also, we depict the solutions for the corre- 
sponding energy densities; see Figure 3. 
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FIG. 1: Numerical BPS solutions to g (r) for the models pTj) 
(dash-dotted red line) and (|27|l (dashed blue line) . The stan- 
dard solution is also depicted (solid black line), for compari- 
son. 

In Figure 1, we present the solutions for the profile 
function g (r) , and we see that both the modified profiles 
reach their vacuum values more slowly than their canon- 
ical counterpart. In this sense, such solutions have a core 
which is greater than the usual Maxwell-Higgs one. Here, 
the conclusion is that, in general, the introduction of a 
non-canonical dynamics allows for the existence of self- 
dual field solutions g (r) with a increased characteristic 
length. Also, we note that the solution for p4| and ((25)) 
goes to its vacuum configuration more slowly than that 
for d^ni) and ((201) • So, beyond the fact of increasing the 
core of the solutions, we note that (pij) increases it more 
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than the ((27)) . We believe that this fact is related with 
the corresponding self-dual Higgs potential, which is of 
sixth-order for the model (|2T|) . and of fourth-order for 
the ([77]) one; see equations ([^ and ([TU|) above. 




r 



FIG. 2: Numerical BPS solutions to a(r). Conventions as in 
FIG. 1. 

In Fig. 2, we depict the numerical solutions for the 
function a (r). Here, as in the Fig. 1, we note that the 
modified solutions go to their vacuum states more slowly 
than the standard Maxwell-Higgs profile, and so such so- 
lutions have a increased characteristic length. This be- 
haviour reinforces the previous conclusion, which states 
that a non-standard dynamics leads to a self-dual field 
solutions with a increased core. We also see that, be- 
yond the fact of increasing the core of a (r) , the model 
pT|) enlarges it more than the (|27|) . Finally, we point 
out the existence of a proeminent plateau in the profile 
related to ([2T|) . near the origin: such structure is also 
present in the self-dual Chern-Simons-Higgs case, which 
is governed by the potential (|33l) . In this context, it is 
interesting to note that the existence of such proeminent 
plateau seems to be closely related to the vacuum mani- 
fold of a sixth-order symmetry breaking potential, since 
it is also present in the modified model (|2T]) . which is 
governed by ([2^ . 

The Fig. 3 encloses the solutions for the energy densi- 
ties of the new BPS states; see equations (^51 and (1^ . 
We see that the profile for the energy density correspond- 
ing to the modified model (|2T|) is quite different from that 
of the standard Maxwell-Higgs theory: in the canonical 
case, the energy density reaches its maximum value in 
the limit r — ^ 0, and it is monotonically decrescent for 
all r. On the other hand, in the modified case (HH, the 
energy density reaches its maximum value at some finite 
distance R from the origin, and it is not monotonically 



decrescent for all values of the independent variable. In 
this context, this difference reinforces our previous con- 
clusion, since such modified behaviour mimics the one of 
the energy density of the usual self-dual Chern-Simons 
theory. 



2,5 
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FIG. 3: Plots of the energy density e(r). Conventions as in 
FIG. 1. 

The solution for the BPS energy density of the modi- 
fied model (P7)) is also depicted in Figure 3, and we see 
that its behaviour is qualitatively similar to that of the 
Maxwell-Higgs model: it reaches its maximum value as 
r goes to 0, and it is monotonically decrescent for all r. 
Here, however, there are two important differences: the 
first one is on the maximum value of the modified profile 
itself, which is smaller than the usual one, and the second 
one is on the characteristic length of the non-standard so- 
lution, which is greater than its canonical counterpart. 

An important issue concerning the study of topological 
structures is the computation of their topological charges, 
which must be conserved. In the present case, i.e., for 
electrically non-charged field solutions of the form © and 
([7]) , this charge is related to the magnetic field generated 
by such solutions themselves. To investigate this issue, 
we introduce the topological current 

= e^^'^^d^Ax , (37) 

which is clearly conserved. Here, e*^^^ = —1. The 0*^- 
component of such current, that is, the topological charge 
density, can be written as 

J° = d^^Ay - dyA^ - B , (38) 

and we see that it is directly proportional to the magnetic 
field B. In this case, the topological charge is given by 

Qt = j BSr = ^B ■ (39) 
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Here, <I>b stands for the flux of the magnetic field. 

According ([7]), ([T7]) and (IT^ . the magnetic flux $b can 
be rewritten in the form 



27rn 

e 



(40) 



and we note that the topological charge Qt is conserved, 
and it is quantized according to the winding number n. 

In this sense, to better specify the field configurations 
studied in this work, a rather natural way is to consider 
their corresponding magnetic fields. In the present con- 
text, this field is given by 



B 



1 da 
r dr 



(41) 



The modified numerical profiles for the magnetic field 
(|¥T|) are plotted in Figure 4. Also for such field, the non- 
usual dynamics introduced earlier leads to BPS solutions 
with a increased core, since the modified profiles have a 
characteristic length which is greater than the standard 
one. In this for the energy densities previously 

depicted, the magnetic field associated to (PT|) is quite 
different from the canonical self-dual one, which reaches 
its maximum value for r — > 0, and is monotonically de- 
crescent ever: the non-standard solution reaches its max- 
imum value at some finite distance from the origin, and it 
is not monotonically decrescent for all r. Such behaviour 
is similar to that of the self-dual magnetic field related 
to the canonical Chern-Simons-Higgs model. 




FIG. 4: Plots of the magnetic field B (r). Conventions as in 
FIG. 1. 

To end this Section, we discuss the solution for the 
magnetic field (jll]) related to the model ([27]). This solu- 
tion is also depicted in Fig. 4, from which we note that 
the behaviour of such magnetic field is quite similar to 



that of the usual Maxwell-Higgs case, since it reaches its 
maximum value as r — >■ 0, and it is monotonically de- 
crescent ever. Here, we point out the behaviour of the 
modified solution, which assures the conservation of the 
topological charge; see eas. (|39)) and pO)) . 



V. ENDING COMMENTS 

In the present letter, we have considered the existence 
of rotationally symmetric BPS solutions in a (2 -I- 1)- 
dimensional space-time. We have investigated such so- 
lutions in a modified self-dual Maxwell-Higgs model en- 
dowed by a non-standard dynamics. Here, the modifica- 
tion was introduced in terms of two non-trivial functions, 
G(|0|) and z/;(|(/)|), which must be positive, in order to 
avoid problems with the energy density (jl5l) of the model. 
So, while G (101) couples with the Maxwell term and, as a 
consequence, changes the dynamics of the gauge field in a 
non- usual way, w (li/)!) couples with the squared covariant 
derivative of the non-charged scalar field, then leading to 
a non-standard dynamics to such field. In this context, 
consistent first order equations were obtained since the 
non-trivial functions G(|(/)|) and ?i;(|(/)|) are constrained 
to the symmetry breaking Higgs potential (|0|) of the 
modified model; see Eq. ([TC)) . 

We have integrated the BPS equations by means of the 
relaxation method, and the numerical results we found 
are depicted in Figs. 1 and 2, for some interesting choices 
to the functions G (|0|) and w (|0|); see eas. (|2T|) and (|27| . 
According these solutions, we conclude that both the 
profile functions g (r) and a (r) are quite sensible to the 
choices made for G(|0|) and lud^l). In particular, it is 
important to reinforce that each of the models (|21l) and 
(|27p is related to a very specific symmetry breaking po- 
tential F (|(/)|) and, as a consequence, such models present 
distinct vacuum manifolds; for details, see eas.([TU| and 

dm. 

Also, using the previous results we found for g (r) and 
a (r) , we have integrated both the energy density e (r) 
and the magnetic field B (r) for the modified models (PT|) 
and ((27)) . and these solutions are depicted in Figs. 3 
and 4, respectively; see eqs.(1221), (UHl) and (gl]). The nu- 
merical analysis reveals that the solutions corresponding 
the model (PT|) behave as those predicted by the Chern- 
Simons-Higgs theory, which reach their maximum val- 
ues at some finite distance from the origin, and are not 
monotonically decrescent for all r. On the other hand, 
the solutions corresponding the model (l?7l) behave as the 
Maxwell-Higgs ones, since they reach their maximum val- 
ues as r — >■ 0, and are monotonically decrescent for all val- 
ues of the independent variable. So, as a consequence, we 
conclude that also the energy density e (r) and the mag- 
netic field B (r) are both sensible to the choices made for 
cm and «; (101). 

We hope that this work may stimulate subsequent 
analysis in the field, concerning mainly the features that 
the modified solutions engender. Also, we point out that 
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the variation on the characteristic length of the modified 
solutions presented here is closely related to the efi'ects 
of anisotropy, which is a feature tipically present in the 
effective field models used to describe the behaviour en- 
gendered by low-energy condensed matter systems. In 
this sense, we point out that such effective models are 
usually based on the Lorentz-breaking ideia, since it in- 
troduces the issue of anisotropy explicity. In the context 
of Lorentz-violating models, such effects where already 
studied by some of us; see, for instance, Ref. [20j . 

A very interesting issue concerns the use of modified 



models to mimic the very same solutions engendered by 
the standard Maxwell-Higgs theory. In particular, such 
issue is now under consideration, and we hope to report 
new results in a near future. 
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